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Abstract 
> 

'sj" ■ Three-jet production in the reaction ep — > eXp has been studied with the ZEUS detector 

^? I at HERA using an integrated luminosity of 42.74 pb~^. The data were measured in 

^ ; the kinematic region 5 < Q^ < iqo GeV^, 200 < W" < 250 GeV and 23 < Mx < 40 

^D ' GeV. The diffractive signal was selected by requiring a large rapidity gap in the outgoing 

proton direction. Jets were reconstructed in the centre-of-mass system of X using the 

^ exclusive /cy-algorithm. A sample of three-jet events in diffraction has been identified. 

g_j ■ Differential cross sections were measured as a function of the jet pseudorapidity and jet 

JL . transverse momentum with respect to the virtual photon-Pomeron axis. The jets going 

g_) I in the Pomeron direction are broader than those going in the virtual-photon direction. 

>-^ ■ This is consistent with models predicting that gluons are predominantly produced in the 

Pomeron direction and quarks in the virtual-photon direction. 
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1 Introduction 

One of the most important results from the ep colhder HERA is the observation that about 
10% of deep inelastic scattering (DIS) events exhibit a large rapidity gap (LRG) between the 
direction of the proton beam and that of the first significant energy deposition in the detector [jl[ . 
These LRG events, ep -^ eXN, result predominantly from the diffractive dissociation of the 
virtual photon, 7*p -^ Xp, and have been interpreted in terms of the exchange of a colour- 
singlet object known as the Pomeron {IP), introduced to describe the energy dependence of 
total cross sections in hadron-hadron scattering 0. 

The HERA data 0, ^ have been analysed in terms of a diffractive structure function, F,p [^, 
defined in analogy with the proton structure function, F2. The results are consistent with a 
"resolved Pomeron" model where F:^ factorises into a Regge-inspired Pomeron flux and a 
Pomeron structure function. The latter has been analysed using the DGLAP evolution 
equations. The extracted parton densities are dominated by gluons. 

Alternatively, a number of QCD models has been proposed to describe diffraction in DIS. 
In these models, the virtual photon dissociates into a qq or qqg state that interacts with the 
proton by the exchange of a gluon ladder. The diffractive cross section can then be formulated 
in terms of contributions from qq and qqg final states. The qqg component is expected to 
dominate the diffractive cross section at high masses, Mx 0, |, |]. 

An important feature of the three-parton final state in both classes of models is that, in 
the rest frame of the system X (7*iP centre-of-mass system for the resolved Pomeron model), 
the qq system tends to populate the virtual-photon hemisphere, while the gluon is emitted in 
the opposite direction (the Pomeron direction in the resolved Pomeron model). In the resolved 
Pomeron model, the gluon is the remnant of the Pomeron left after the boson-gluon fusion 
process has produced a qq pair; as a remnant, it follows the original Pomeron direction. In 
dipole models, to leading order, the gluon in the qqg fluctuation from the virtual photon either 
couples directly to the exchanged gluon ladder or emerges from it, depending on the reference 
frame. In this case also, therefore, it follows the direction of the gluon ladder (= Pomeron in 
resolved Pomeron model) and the qq pair carries on in the direction of the virtual photon. 

Studies of the hadronic final state and jet production in diffractive interactions have been 
presented by both HI |]lOl and ZEUS ]ll|. In this paper, the jet-like structure of the hadronic 



final state is studied with the aim of isolating three-parton final states by selecting three-jet 



events. As was already observed in e~^e~ experiments |jT2[, three-jet topologies can be clearly 
separated only for centre-of-mass energies above about 20 GeV. The measurements are thus 
performed in a restricted phase-space region populated by high-mass diffractive systems. Such 
systems are a small fraction of the total diffractive cross section. They are separated from a 
substantial background of non- diffractive processes by requiring a large rapidity gap. A search 
is made for jets using an algorithm that allows the study of jet configurations aligned with 
respect to the 7*iP axis and ensures good parton-hadron correlations. The properties of the 
three-jet final states are used to test the expectation that the gluon from the dominant qqg 
contribution is emitted in the Pomeron direction. Cross sections are presented as a function of 
the jet pseudorapidity and jet transverse momentum with respect to the ■j*]P axis. In addition, 
the internal structure of the jets is measured for jets found in either the photon or the Pomeron 
directions, with the aim of distinguishing between quark- and gluon-initiated jets. The different 
measurements are compared to models of diffraction in DIS. 



2 Experimental setup 

The data were collected with the ZEUS detector at HERA during 1998-2000 when HERA 
collided 27.6 GeV electrons or positronsQ on 920 GeV protons. A total of 7.96 ± 0.14 pb~^ of 
e~p data and 34.78 ± 0.87 pb~^ of e+p data was used. A detailed description of the ZEUS 



detector can be found elsewhere ||T3|, |1J]. A brief outline of the most relevant components for 
this analysis is given below. 

The ZEUS high-resolution uranium-scintillator calorimeter (CAL) [|15| covers 99.7% of the 
total solid angle. It is divided into three parts with respect to the polar angleQ as viewed from 
the nominal interaction point: forward (FCAL), rear (RCAL) and barrel (BCAL) calorimeters. 
Each part is subdivided longitudinally into electromagnetic (EMC) and hadronic (HAC) sec- 
tions. The CAL energy resolution, as measured in test-beams, is o{E)/ E = 0.18/ Je (GeV) 



for positrons and a{E)/E = 0.35/ JE (GeV) for hadrons. 

Charged particles were tracked by the central tracking detector (CTD) |l^, which operates 
in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD is a drift 
chamber consisting of 72 cylindrical layers, organised in 9 superlayers, covering the region 
15° < 9 < 164°. The transverse momentum resolution for full-length tracks is cr{pt)/pt = 
O.OObSpt © 0.0065 © 0.0014/pt, with pt in GeV. 

In 1998, a Forward Plug Calorimeter (FPC) |]l^ was installed in the 20 x 20 cm^ beam hole 
of the FCAL with a small hole of radius 3.15 cm in the centre of the FPC to accommodate the 
beam pipe. This increased the forward calorimetric coverage by about 1 unit of pseudorapidity 
to —3.8 < ?7 < 5. The FPC consisted of a lead-scintillator sandwich calorimeter divided into 
electromagnetic and hadronic sections which were read out separately by wavelength-shifting 
fibres and photomultipliers. The EMC (HAC) sections were segmented into cells of 2.4 x 2.4 cm^ 
(4.8 X 4.8 cm^) area and constituted a total of 76 read-out channels. The energy resolution was 
(Je/E = 0.41/vE' © 0.062 and (Je/E = 0.65/vE' © 0.06 for positrons and pions, respectively, 
where the energies are measured in units of GeV |[T7| . 

The small-angle rear tracking detector (SRTD) [|l8l was used to measure the position of 
those positrons scattered at a sufficiently small angle to strike it. It consists of two planes 
of scintillator strips attached to the front face of the RCAL. The SRTD signals resolve single 
minimum-ionizing particles and provide a position resolution of 0.3 cm. The position of those 
positrons falling outside the SRTD acceptance was measured in the CAL. 

The luminosity was measured via the positron-proton bremsstrahlung process, ep -^ e'jp, 
using a lead-scintillator calorimeter (LUMI) |jl9[ located at Z = —107 m in the HERA tunnel. 



^Hereafter "positron" is used to refer to both electron and positron beams. For the Q^ range studied, e~p 
and e'^p scattering are assumed to give identical results since contributions from Z'^ exchange are negligible. 

^The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the proton 
beam direction, referred to as the "forward direction" , and the X axis pointing left towards the centre of HERA. 
The coordinate origin is at the nominal interaction point. The pseudorapidity is defined as 77 = — ln(tan|), 
where the polar angle, 9, is measured with respect to the proton beam direction. 



3 Kinematics 

The kinematics of inclusive deep inelastic scattering e{k) p{P) — > e{k') + anything is described 
by any two of the following variables: 

Q = -q = -[k~k) , x= —- — , y = — — , W = + m , 

2F ■ q F ■ k X ^ 

where k and k' are, respectively, the four-momenta of the initial and final positrons, P is the 
initial proton four-momentum, y is the fraction of the energy transferred to the proton in its 
rest frame, mp is the proton mass and W is the 7*p centre-of-mass energy. 

For the description of the diffractive process ep -^ eXp, in addition to the invariant mass of 
the system X, Mx-, two further variables are introduced: 






X 



In resolved Pomeron models, where the interaction is described as the exchange of a particle- 
like Pomeron, xp is the fraction of the proton momentum carried by the Pomeron and /3 is the 
fraction of its momentum carried by the parton within it that is probed by the virtual photon. 

4 Monte Carlo simulations 

Three Monte Carlo (MC) models for the diffractive deep inelastic process ep -^ eXp have been 
considered in order to describe the diffractive hadronic final state. 

The RAPGAP 2.08/06 MC generator |^ implements the Ingelman-Schlein factorisable 
model for the Pomeron. The matrix elements for the 0{a) process (eg —>■ eg) and for the 
0{aas) processes (eg — > eqg, eg -^ eqq) are included. To avoid divergences in the matrix 
elements for the 0{aas) processes for massless quarks, a cut of pf, > 3 GeV^ was applied, 
where Pt is the transverse momentum of any of the outgoing partons with respect to the 
photon direction in the centre-of-mass frame of the hard scattering. The relative contribution 
of each component is determined by the quark and gluon densities within the Pomeron, which 
evolve according to the DGLAP equations. The HI parameterisation [Q for the Pomeron 
parton densities with a 'hard gluon' carrying > 80% of the momentum at Qo = 3 GeV^ was 
used. Two different RAPGAP MC samples were generated, in which the higher-order QCD 
corrections were treated in two different ways: either using the colour-dipole model (CDM) |2l|] 
as implemented in ARIADNE 4.08 [22| or using parton showers (PS) as implemented in LEPTO 
6.1 p3[. The Lund string-fragmentation scheme ||2^ as implemented in JETSET 7.4 [^ was 
used for hadronisation. First-order electroweak corrections are taken into account via the 



HERACLES 4.6 program 12| 



The RIDI 2.0 MC generator ^\ implements the diffractive dissociation of the virtual pho- 
ton into gg and qqg final states following the approach of Ryskin ^ , in the framework of the 
leading logarithmic approximation of pQCD. Contributions from transversely and longitudi- 
nally polarised photons are included. The predicted cross section is proportional to the square 
of the gluon momentum density in the proton. The CTEQ4M parameterisation [^ for the 



parton densities was used. Final-state parton showers and fragmentation were implemented 
using JETSET 7.4. First-order radiative corrections were also included. 



The SATRAP MC generator [^ is based on the Golec-Biernat and Wiisthoff model |^ 
In this model, the diffractive cross section is expressed as the convolution of the photon wave- 
function, calculated for qq and qqg colour dipoles and including the contributions of both 
transversely and longitudinally polarised photons, with the dipole cross section for scattering 
on the proton. The parameters of the model were tuned to describe the total 7*p cross section as 
measured at HERA. Final-state parton showers and fragmentation were simulated via JETSET 
7.4. First-order radiative corrections were taken into account via the HERACLES 4.6 program. 

During the course of this analysis, it became clear that neither RIDI nor SATRAP gave a 
reasonable description of the data. As discussed in Section 8, this was traced to inadequacies 
in the implementation of the modelling of higher-order QCD processes in these models. Most 
notably, no initial-state parton cascades were included, and the final-state QCD radiation from 
the gluon in the dominant qqg contribution was suppressed. A new implementation of higher- 
order QCD processes in SATRAP was carried out |3l|, ^, in which the CDM approach was 
implemented in a similar fashion to that in RAPGAP. This model is referred to as SATRAP- 
CDM. 

Processes where the proton dissociates into a system N , ep — *■ eXN, were generated using 



the EPSOFT 2.0 MC generator [Q. This simulates events based on the triple- Regge formalism 
Q] and an Mpf distribution as measured for diffractive dissociation in pp scattering. Samples 
with different parameters for the simulation of the system N were considered. 

Non- diffractive DIS background was generated including first-order QED radiative correc- 
tions using LEPTO 6.5 interfaced to HERACLES 4.5 via DJANGOH 1.1 [§§. The CTEQ4M 
set of proton parton densities was used. The colour-dipole model was used to simulate parton 
cascades and fragmentation was performed using JETSET 7.4. 

The RAPGAP-CDM model provides a good description of most of the measured distribu- 
tions and was used to study the accuracy of the kinematic variables and jet reconstruction, 
the efficiency for selecting events and the corrections for detector and resolution effects. The 



detector simulation was based on the GEANT 3.13 program ||3^ 



5 Event reconstruction and selection 

The variables Q^ and W were determined from the information of the scattered positron. 
The hadronic final-state system, X, resulting from the dissociation of the virtual photon, was 



reconstructed using "energy fiow objects", hereafter denoted as EFOs fS^], assumed to have the 
mass of the pion. These EFOs combine information from charged tracks, as measured in the 
CTD, and energy clusters measured in the CAL and FPC. In the kinematic region considered 
in this analysis, the invariant mass Mx was reconstructed on average with a systematic shift of 
—14% and a relative resolution of 10%, while the hadronic three- vector was reconstructed with 
a resolution of about 1 GeV for the X and Y components and ~ 2 GeV for the Z component. 
The four-momentum of the virtual photon, q, was determined with no significant bias and with 
a relative resolution of about 10% for the X and Y components and ~ 7% for the Z and energy 
components. 

The triggering and online event selections were similar to those used for the ZEUS measure- 



ment of the structure function F2 |33]- The offline selection criteria, applied to select diffractive 
DIS events, required: 



a scattered positron candidate (identified |3^ via its pattern of energy deposition in the 



calorimeter) with energy, E'^, greater than 10 GeV and y < 0.95. For an incoming positron 



of energy E^,, y was estimated from the energy and the polar angle, 6'^, of the scattered 
positron by: 

The impact position of the scattered positron in the RCAL was required to be outside a 
box of 26 X 20 cm^ around the beam pipe to ensure a fully contained positron candidate. 
In addition, candidates within the CTD acceptance region, but without a matched charged 
track reconstructed in the CTD, were rejected; 

• 40 < (5 < 70 GeV, where S = J2i Ei{l — cos6'i) and the sum runs over all calorimeter cells, 
i.e. including those belonging to the scattered positron. This cut, together with the previ- 
ous requirements on the positron candidate, eliminated background from photoproduction 
events and beam-gas interactions producing fake scattered positrons; 

• a reconstructed vertex with \Z\ < 50 cm and at least two associated tracks; 

• the presence of a large rapidity gap in the outgoing proton direction, to select diffractive 
events. Events with /y™**^ < 3.0 were selected, where r^™^^ is the pseudorapidity of the 
most-forward cluster with energy above 400 MeV, as measured in the CAL and FPC; 

• at least three reconstructed EFOs in the hadronic final state to ensure that the event 
plane could be determined. 

The measurements were limited to the following kinematic range: 

• 5 < g2 < 100 GeV^; 

• 200<W < 250 GeV; 

• 23 < Mx < 40 GeV, where Mx was corrected for energy losses (see Section 6.1); 

• xjp < 0.025. 

In this region, the diffractive hadronic final state X is well contained in the central detector, 
thus minimising energy losses through the forward and rear beam-pipe holes. A good separa- 
tion of the diffractive signal from the non-diffractive background was achieved: after the r/™*^^ 
requirement, the background was typically around 3% and at most 10%, as estimated from MC 
studies. After all cuts, 7175 events remained. 



6 Three-jet search 

A search for jets was performed in the centre-of-mass system of the observed hadronic final 
state, X, using the exclusive A;r-algorithm [^. All EFOs in the event were used as input to 
the jet-finding procedure. For each pair of particles (EFOs) i and j in an event, the quantity 

2-min(E2,E|)(l-cos%) 



Vij 



Ml 



was computed, where Ei is the energy of the i*'^ particle and 9ij is the angle between particles 
i and j. The pair of particles with the smallest value of yij was replaced by a pseudo-particle 



or cluster. The four-momentum of the cluster was determined using the "E-recombination 
scheme" : 

^cluster = Ei + Ej , Pfc^'"^*^-- = Pi + Pl {k = X, Y, Z) . 

This clustering procedure was repeated until all yij values exceeded a given threshold, ycut, and 
all the remaining clusters were then labelled as jets. For this recombination scheme, the jets 
become massive and the total invariant mass, Mx, coincides with the invariant mass of the 
three-jet system. As applied in the centre-of-mass system, this algorithm produces at least two 
jets in every event. The same jet search procedure was applied to the final-state hadrons for 
simulated events. 

Figure |T] shows the measured rijet fractions for rijet = 2,3 and > 3 as a function of the jet res- 
olution parameter, ycut, in the region 0.01 < ^cut < 0.1. The rate of three-jet production varies 
from 47% at ?/cut = 0.01 to 8% at T/cut = 0.1. The measured jet fractions were compared with 
diffractive MC models. All the models show trends similar to the data but only SATRAP-CDM 
provides a good description of the measurements. The predicted rijet fractions are sensitive to 
the particular implementation of higher-order QCD corrections and hadronisation processes in 
each MC model. The differences observed between the predicted jet fractions from RAPGAP- 
CDM and RAPGAP-PS for ?/cut ~ 0.03 indicate where the dependence on the model used for 
the simulation of higher-order QCD contributions starts to be sizeable. 

Good parton-hadron correlation in three-jet production was observed in the MC simulation 
for Mx > 20 GeV. In this mass range, values for ?/cut around 0.05 gave the smallest hadronisation 
corrections. The parton-hadron jet correlation was studied using a MC sample with three jets 
found in the final state at both parton and hadron levels. The angular distance between 
each parton and the corresponding jet of hadrons was such that the partons were always well 
contained within the size of the hadron jets, which have a typical angular radius larger than 
0.5 radians. 

6.1 Jet reconstruction and selection 

The reconstruction of the jet variables was studied using MC events with three jets at both 
hadron and detector levels. Matched pairs of jets were selected and used to study the jet 
reconstruction in the detector. The characteristics of the reconstructed jet quantities were 
found to be as follows: 

• the polar angle of the jet, 6'*,^, defined with respect to the 7*iP axi^, was reconstructed 
with no significant systematic bias and a resolution of ~ 6° and ~ 8° for the two most 
energetic jets and the least energetic jet in the event, respectively; 

• the measured invariant mass, Mx, and jet energies, E*^^, were corrected for the energy 
loss by a common factor of 1.14, determined from MC studies. In the kinematic range 
considered, the jet energy was reconstructed with a typical resolution of 10% and the jet 
energy correction was independent of the jet direction; 

• the transverse momentum of the jets, -Rfjef measured with respect to the 7*iP axis, was 
reconstructed without significant bias and with a resolution of 0.7 GeV; 



^In the centre-of-mass system, the Pomeron defines the positive Z direction. 



• the jet pseudorapidity, rjl^, measured in the laboratory frame (see below), was recon- 
structed without significant bias and with a resolution of 0.1 to 0.2 units. 

In addition, it was checked that, in the Mx region considered and for a fixed value of 
2/cut = 0.05, the jets, defined in the centre-of-mass frame, always had energy significantly above 
the noise level in the CAL. The jet energies in the laboratory frame, Ell-^, computed by boosting 
the jet four-momenta in the centre-of-mass system back to the laboratory frame, were typically 
above 4 GeV, thus ensuring a well understood jet-energy reconstruction for all jets. 

The sample of three-jet final states was defined by using a fixed ?/cut = 0.05 and the jets 
were required to be in the range |?7i^J^| < 2.3. The variable r]l'^ was computed in the same way 
as -Slab- The effect of the inclusion of jet masses in the determination of r]ll^ was found to be 
negligible. This procedure selected events with jets reconstructed in the central part of the 
detector, which are not affected by the r/™'™ cut applied in the forward region and for which a 
good correlation between the jet before and after detector effects was preserved. After the r^J^j, 
cut, a final sample of 891 events with exactly three identified jets was obtained. 

6.2 Three-jet topology variables 

Neglecting jet mass, the topology of a three-jet final state can be described using the fractional- 
momentum variables 

^ " \^ i 

y. ipJ: 



x^ = ^T^' |pf I = v(^5)^ + (pi^y + (pSr > (1) 



where P]^* (i = 1 — 3) denotes the i^^ jet three-momentum, determined in the centre-of-mass 
frame, and by construction X1+X2 + X3 = 2. The jets were sorted according to their momenta 
in such a way that 

Xl> X2 > X3 . 

Two of the Xi variables in Eq. (|l]) are sufficient to describe the jet topology in the final state. 
Following the studies of three-jet production in e~^e~ collisions at PETRA [EO, O, these quan- 



tities were chosen to be Xi and $, = {x2 — x^)/xi, as suggested by Ellis and Karliner [^ 



6.3 Jet-shape reconstruction 

The internal structure of the jets was studied in terms of the differential jet shape, defined in 
the centre-of-mass system as the average of the fraction of the jet energy which lies inside an 
annulus of inner angular distance (p — 5ip/2 and outer angular distance ^p + 5ip/2 from the jet 
axis. EFOs were used to reconstruct the jet shape defined as 

EFOs, 1 ^ae;,(^±^^/2) 



p''°^(v') = 7:7^E 



^^ ^jets ^ E^ 



jet 



where A£'jg(.((^ ± Sip/ 2) denotes the sum of the energies of the EFOs belonging to a given jet 
whose angular distance to the jet axis is within the range {p — 5p:>/2, p} + 5p/2) with 5p = 0.2 
radians; A'jets is the number of jets. For MC events, the same jet-shape definition as used above 
for the EFOs was applied to the final-state hadrons. 



7 Unfolding and systematic studies 

The measured distributions were corrected bin-by-bin for acceptance and migrations using 
RAPGAP-CDM. The measurement was carried out in the region 5 < Q^ < 100 GeV^, 
200 < W < 250 GeV, 23 < Mx < 40 GeV and xp < 0.025, for diffractive events with 
^hadron < ^.O and | ^^j^b | hadron < 2.3, whcrc ^h^d^n ^^ the pscudorapidity of the most-forward 
final-state hadron with energy above 400 MeV. The measurements presented here thus cover 
a very restricted phase-space region. According to RAPGAP, in the kinematic region con- 
sidered, less than 10% of all diffractive events have //hadron < 3.0 and three jets in the final 
state with |?7iablhadron < 2.3. The fraction of three-jet events selected by the //h^dmn < 3.0 and 
l^iablhadron < 2.3 requirements was estimated using RAPGAP to be of the order of 30%. 

A detailed study of the main sources contributing to the systematic uncertainties of the 
measurements was performed |31]]. The different contributions are listed below, with a typical 
value quoted for each: 

• the uncertainty in the hadronic final-state energy scale was taken into account by varying 
the GAL (FPG) energy by ±3% (±10%) in the MG, after removing the energy deposits 
belonging to the scattered positron. This introduced an uncertainty in the measured 
distributions of about 5%; 

• the selected window in the reconstructed Mx was varied by ±10%, both in data and sim- 
ulated events, according to the estimated resolution in the reconstruction of this quantity. 
The measured cross sections varied by 5%; 

• the uncertainty in the positron energy scale was included by varying the positron energy 
by ±2% in the MG. The effect on the r]*,^ and -Pyjet distributions (see Section 8) was 
typically below 3%; 

• the box cut applied to the position of the scattered positron in the RGAL was varied by 
±5 mm (both in data and simulated events) to account for the uncertainty in simulating 
the measured positron position. The effect was typically below 1%; 

• the selected range in the reconstructed jet pseudorapidity in the laboratory frame, |?7i^b| < 
2.3, (see Section 6.1) was varied by 0.2 units around the nominal value, both in data and 
MG simulation, according to the estimated resolution in the reconstruction of the jet 
direction. This introduced an uncertainty on the measured cross sections of about 2%; 

• the measured distributions were corrected to the hadron level using RAPGAP-PS instead 
of RAPGAP-GDM to account for the uncertainty in the simulation of higher-order QGD 
processes in the MG model. The effect varied between 5% and 10% and was the dominant 
systematic uncertainty. 

The final systematic uncertainty was computed by adding the different contributions in 
quadrature. The contribution from proton-dissociative processes, ep — > eXN, with Mx < 
3 GeV, where the system A^ escaped through the forward beam hole in the FPG, was estimated 
using EPSOFT to be (16 ± 5)% and was subtracted. The overall normalisation uncertainty of 
5.4% in the measured cross sections, coming from the uncertainty on the proton-dissociative 
contribution and the 2% error in the luminosity determination, is not included in the figures. 
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8 Results 

Multi-jet production in the centre-of-mass system of the hadronic final state, X, is presented 
in the kinematic region defined in Section 7 using the exclusive /cT-algorithm with y^ut = 0.05. 
The measurement^ were corrected for electroweak radiative effects to the Born level. 

8.1 Three-jet topology 

The topology of the three-jet final state was studied using Xi, the fractional-momentum vari- 
ables (see Section 6.2). Figure ||(a) shows the distribution of the three-jet sample in the (^, Xi) 
plane. As indicated in the figure, different regions in the plane correspond to different three-jet 
topologies. All configurations are present, including those for which the three jets have similar 
energies. The distribution of the events shows that, for the given ?/cut, the favoured three-jet 
topologies are those for which one leading jet carries approximately twice the energy of either 
of the other two jets (^ ~ 0.1, Xi ~ 0.9). The normalised differential cross sections for three-jet 
production as a function of Xi and ^ are presented in Figs. 0(b) and 0(c), compared to the MC 
predictions. The measured Xi distribution has a peak structure with a maximum at Xi ~ 0.9, 
which is determined by the ?/cut value used in defining jets since the position of the peak ap- 
proximately corresponds to 1 — 2?/cut- The shapes of the measured cross sections are reproduced 
by all three MC models. This indicates that they provide a good description of the dynamics 
of the three-jet emission, which dictates the population of the (^, xi) plane. 

The energy fiow, as measured in the event plane defined by the two most-energetic jets in 
the event, was studied as a function of the azimuthal angle, 0*. The origin of 0* was defined 
to be along the direction of the most-energetic jet with positive 0* increasing in the direction 
of the second most-energetic jet. The energy fiow was measured by projecting the particle mo- 
menta onto the event plane and summing, in each 0* bin, the particle energies. The sums were 
normalised by dividing by Mx- In Fig. ^ the normalised energy fiow is presented for different 
regions of the (^, xi) plane. A clear three-jet structure is observed in all regions. The measure- 
ments are compared with the RAPGAP and SATRAP-CDM predictions. RAPGAP-CDM, 
RAPGAP-PS (not shown but similar to RAPGAP-CDM) and SATRAP-CDM all provide a 
very good description of the measured fiows for all (^, Xi) regions and 0* ranges. 

8.2 Cross section as a function of 77 ^^^ 

Figure ^(a,b) shows the differential cross section as a function of the jet pseudorapidity, da/dri*^^, 
measured with respect to the 7*iP axis. All jets in each event were included. The total cross 
section for three-jet production, obtained by integrating the differential measurement, is: 

""'''*' ^\' I Wt'^'^^"'^ ^^'^ ^ 0.5(stat.)li:?(syst.) pb. 

For this measurement, the overall normalisation uncertainty of 5.4% (see Section 7) is com- 
bined in quadrature in the quoted total systematic uncertainty. The measured cross section 
corresponds to 32.2% ± 0.9%(stat.)ti;6%(syst.) of the diffractive cross section for two or more 

■*The measurements are tabulated in Tables 1-7. 



jets. This rate is consistent with being independent of Q^ and Mx- 

Figure ^(a) shows the measured differential cross section compared to the RAPGAP-CDM 
and RAPGAP-PS predictions. Both implementations of RAPGAP provide a good description 
of the shape of the cross section but underestimate it by about 20%. The observed deficit in 
the absolute normalisation predicted by RAPGAP may be attributed to the absence of next-to- 
leading-order terms in the computed matrix elements, which are only approximated via parton 
cascades. 

In Fig. ^(b) the measured differential cross section is compared to the SATRAP, SATRAP- 
CDM and RIDI predictions. For the latter, the contribution from the hadron jet associated with 
the gluon from the qqg final state is shown separately. The SATRAP and RIDI predictions are 
found to be more than a factor of two below the measured cross sections and show a double-peak 
structure not observed in the data. This tendency is particularly pronounced in RIDI. This 
reflects the presence of dominant topologies in these models in which the jet in the Pomeron 
direction is emitted almost coUinear to the 7*iP axis and for which the second and third jets are 
found in the photon hemisphere. The observed shape of the SATRAP and RIDI predictions 
is a consequence of the inadequacies of the implementation of the modelling of higher-order 
QCD processes in these MC models as discussed in Section 4. SATRAP-CDM describes the 
shape of the measured distribution but also underestimates the normalisation by about 20%, 
which, as for RAPGAP, may be attributed to the absence of next-to-leading-order and some 
leading-order terms in the computed matrix elements. 

8.3 P^ ■ . distribution for the most-forward jet 

The differential cross section as a function of the jet transverse momentum, -Pfjet) measured 
with respect to the 7*iP axis for the most-forward jet is presented in Fig. ^(c,d). The mea- 
surement was performed in the region P-^jet > 1 GeV, where a good parton jet-hadron jet 
correlation is preserved, and for which the measurement is well defined, given the estimated 
resolution in the reconstruction of -P^jet i^ the detector (see Section 6). The measured cross 
section is dominated by configurations with P^ ■ ^ below 3 GeV but contains a long tail towards 
larger -Pfjet values, reflecting the presence of contributions from higher-order QCD processes. 

Figure ^(c) shows the measured cross section compared to RAPGAP-CDM and RAPGAP- 
PS. The RAPGAP predictions underestimate the measured cross section and have a slightly 
different shape than the data. RAPGAP-CDM provides a better description of the data than 
RAPGAP-PS. The observed difference between these two predictions indicates the sensitivity 
of the -Rfjet distribution to the approach used for the simulation of higher-order QCD processes 
in the MC. 

Figure ^(d) shows the measured -Pfjet distribution compared to the predictions from SATRAP, 
SATRAP-CDM and RIDI. The P^j^t distribution predicted by SATRAP and RIDI is steeper 
than that of the data, indicating that, in these MC models, the large Pf jet configurations are 
missing. The cross section in these models is dominated by configurations strongly aligned with 
respect to the 7*iP axis in the forward region. As already mentioned, this is most probably a 
consequence of the inadequacies of the implementation of the modelling of higher-order QCD 
processes in these models. In SATRAP-CDM, the shape and normalisation of the P^jet distri- 
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bution is in better agreement with the data. 

In all models, the gluon in the qqg final state is emitted with very small transverse momentum 
with respect to the 7*iP axis. The contributions from higher-order QCD processes are expected 
to be important and to infiuence the kinematics of the corresponding hadron jet which, due to 
string effects, tends to have smaller pseudorapidity (larger transverse momentum) compared to 
the initial parton. Given the |^iab|hadron < 2.3 requirement applied in defining the measured cross 
sections, the absolute predictions from these models depend on an appropriate treatment of 
these higher-order QCD emissions and a proper description of the forward hadron-jet dynamics 
(as shown in Fig. ^(d)). This is particularly clear in the case of SATRAP, which fails to 



reproduce the three-jet data but provides a good description of inclusive diffractive DIS ||29 



8.4 Jet shapes 

The differential jet shapes, p(v5), were measured for the most-forward and most-backward jet 
in three-jet events, where the forward region is defined by the Pomeron direction. The mea- 
surements were performed in the 7*iP centre-of-mass frame for jets with energy, E^^^^, above 
9 GeV. This requirement removes the contributions from the least-energetic jet in the event 
which, in the M^-range considered and for the given i/cut, has energy below this threshold and 
for which hadronisation effects on the jet shape are expected to be important. Figure ^ shows 
the measured jet shapes compared with the RAPGAP-CDM predictions (the SATRAP-CDM 
predictions, not shown, are almost identical). The jet in the Pomeron direction is broader than 
the jet in the photon direction. These measurements are described by models in which a gluon 
populates the Pomeron hemisphere and a quark is found in the photon direction. The energy 
distribution of the two jets is similar, indicating that the observed difference in the measured 
jet shapes can be ascribed to the different nature of the initial partons. Moreover, the difference 
between the jet shapes in the Pomeron and photon directions is qualitatively similar to that be- 
tween quark- and gluon-initiated jets as measured in e^e~ experiments [^]. This measurement 
therefore supports the validity of the picture where the three-body final state is dominated by 
a qqg configuration with the gluon preferentially emitted in the Pomeron direction. 

9 Summary and conclusions 

Multi-jet production has been studied in the kinematic region 5 < Q^ < 100 GeV^, 200 < 
W < 250 GeV, 23 < Mx < 40 GeV and xp < 0.025, for diffractive events with r^S^^ < 3.0 
and Ir/i^blhadron < 2.3. A study of three-jet production, using the exclusive fc^-algorithm with 
1/cut = 0.05 in the centre-of-mass system of the hadronic final state, has been presented. The rate 
of three-jet production has been measured to be 32.2% ± 0.9%(stat.)lJ;6%(syst.) of the diffrac- 
tive jet cross section, and is consistent with being independent of Q^ and Mx- The topologies 
of the three-jet final states have been studied in terms of the jet fractional-momentum variables 
S, and Xi\ all possible topologies have been observed. The distribution of the events in the 
(^, Xi) plane shows that, for a given T/cut, the favoured topologies are those for which one lead- 
ing jet carries approximately twice the energy of either of the other two jets (^ ~ 0.1, Xi ~ 0.9). 
Differential cross sections for three-jet production have been measured as a function of the 
jet pseudorapidity, rj^^^, and jet transverse momentum, -Pfjet) "with respect to the 7*iP axis. 
The total cross section for three-jet production, determined by integrating the ri*^^ distribution, 
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is a^i^^^ = 14.1 ± 0.5(stat.)lj;;5(syst.) pb. The -Rfjet distribution for the most-forward jet is 
dominated by configurations with -Rrjet below 3 GeV but exhibits a long tail towards larger 
values, reflecting the presence of higher-order QCD contributions. The differential jet shapes, 
p((/9), of the two most-energetic jets in three-jet events have been measured. The jet in the 
Pomeron direction is broader and consistent with the shape of gluon-initiated jets as predicted 
by diffractive MC models. 

The data are broadly consistent with models in which the hadronic final-state is dominated 
by a qqg system with the gluon preferentially emitted in the Pomeron direction. Such config- 
urations are predicted both by resolved Pomeron models (as implemented in RAPGAP) with 
a Pomeron dominated by gluons, and by models where the virtual photon dissociates in a qqg 
system which interacts with the proton via the exchange of a gluon ladder (as implemented 
in SATRAP and RIDI). RAPGAP provides a reasonable description of the measured distri- 
butions, although it underestimates the measured cross section by about 20%. This deficit 
may be attributed to the absence of next-to-leading-order corrections, which are only included 
approximately via parton cascades. The predictions from SATRAP and RIDI fail to describe 
the shape and normalisation of the 77^ and P-r \ct distributions. However, it has been shown 
that an improved implementation of the modelling of higher-order QCD processes in SATRAP- 
CDM gives reasonable agreement with the data. Thus, the study of three jets in diffraction 
constitutes a very sensitive method to investigate the dynamics of diffractive deep inelastic 
scattering. 
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Figure 1: Measurement of the njet-fraction, for njet = 2,3 and > 3, as a function of the 
exclusive /cT-jet resolution parameter, |/cut- The MC expectations from RAPGAP-CDM (solid 
lines), RAPGAP-PS (dashed-dotted lines) and SATRAP-CDM (dashed lines) are shown. The 
vertical arrow at ycut = 0.05 indicates the value used in the subsequent analysis. 
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Figure 2: (a) Distribution of the three-jet sample in the (^, Xi) plane. The different topologies 
within the plane are indicated, where the length of the arrow represents the momentum of the 
jet as computed in the centre-of-mass frame. (b,c) Normalised differential cross sections as 
a function of Xi and ^ in three-jet production. The thick error bars indicate the statistical 
uncertainties and the thin error bars indicate the statistical and systematic uncertainties com- 
bined in quadrature. The MC expectations from RAPGAP-CDM (sohd hue), RAPGAP-PS 
(dashed-dotted line) and SATRAP-CDM (dashed line) are shown. 
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Figure 3: Energy-flow distribution in the three-jet plane, normalised to the total invariant 
mass, as a function of the azimuthal angle, (f)*, in the different regions of the (1^, xi) plane. The 
azimuthal angle is defined to run from the first to the second most-energetic jet in each event. 
The thick error bars indicate the statistical uncertainties and the thin error bars indicate the 
statistical and systematic uncertainties combined in quadrature. The MC expectations from 
RAPGAP-CDM (solid lines) and SATRAP-CDM (dashed lines) are shown. The characteristic 
topology of the three-jet configuration is indicated in the top right corner of each plot. 
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Figure 4: Differential cross section for three-jet production as a function of (a,b) the jet 
pseudorapidity, da/drj*^^, of the three jets; (c,d) the jet transverse momentum, -PTjef ^^ the 
most-forward jet. The Pomeron direction defines the forward hemisphere. The thick error 
bars indicate the statistical uncertainties and the thin error bars indicate the statistical and 
systematic uncertainties combined in quadrature. The overall normalisation uncertainty of 
5.4% was not included in the figures. The MC expectations are shown from (a,c) RAPGAP- 
CDM (solid line) and RAPGAP-PS (dashed-dotted line) and (b,d) SATRAP (dashed line), 
RIDI (dotted line) and SATRAP-CDM (solid line). For RIDI, the contribution from the gluon- 
initiated jet is shown separately in (b) (shaded area). 
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Figure 5: The differential jet shape, p{ip), for the most-forward and most-backward jets with 
E*^^ > 9 GeV in three-jet events, where the Pomeron defines the forward direction. The thick 
error bars indicate the statistical uncertainties and the thin error bars indicate the statistical 
and systematic uncertainties combined in quadrature. The MC expectations from RAPGAP- 
CDM are shown. 
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A'jct fractions versus ycut 




iVjct = 2 


iVjct = 3 


A^jct > 3 


ycut 


rate 


(stat. syst.) 


rate 


(stat. syst.) 


rate 


(stat. syst.) 


0.01 


0.232 


+0.017 
-0.010 


0.471 


+0.008 
-0.012 


0.297 


+0.007 
-0.009 


0.02 


0.503 


+0.023 
-0.016 


0.420 


+0.016 
-0.024 


0.077 


+0.004 
-0.004 


0.03 


0.655 


+0.021 
-0.012 


0.320 


+0.011 
-0.022 


0.025 


+0.002 
-0.002 


0.04 


0.750 


+0.012 
-0.010 


0.241 


+0.010 
-0.012 


0.009 


+0.001 
-0.001 


0.05 


0.809 


+0.011 
-0.009 


0.187 


+0.009 
-0.011 


0.004 


+0.001 
-0.001 


0.06 


0.850 


+0.011 
-0.007 


0.149 


+0.007 
-0.011 


0.001 


+0.001 
-0.000 


0.07 


0.878 


+0.009 
-0.006 


0.122 


+0.006 
-0.008 


- 


- 


0.08 


0.902 


+0.007 
-0.006 


0.097 


+0.005 
-0.007 


- 


- 


0.09 


0.919 


+0.006 
-0.004 


0.081 


+0.004 
-0.006 


- 


- 


0.10 


0.934 


+0.004 
-0.004 


0.066 


+0.004 
-0.004 


- 


- 



Table 1: Measurement of the njct-fraction, for rijct 
kx jet resolution parameter, |/cut- 



2, 3 and > 3, as a function of the exclusive 



1 da 
a dxi 


dxi 


1 da 
a dx\ 


(stat.) 


(stat. syst.) 


[0.650,0.720] 


0.288 


±0.056 


+0.098 
-0.091 


[0.720,0.790] 


1.639 


±0.162 


+0.226 
-0.345 


[0.790,0.860] 


4.869 


±0.277 


+0.387 
-0.656 


[0.860,0.930] 


6.339 


±0.319 


+0.668 
-0.428 


[0.930,1.000] 


1.150 


±0.150 


+0.319 
-0.266 



Table 2: Normalised differential cross section as a function of Xi. 



\da 
adi 


di 


1 da 
a d£, 


(stat.) 


(stat. syst.) 


[0.0,0.2] 


2.189 


±0.106 


+0.196 
-0.169 


[0.2,0.4] 


1.516 


±0.093 


+0.179 
-0.123 


[0.4,0.6] 


1.153 


±0.086 


+0.114 
-0.225 


[0.6,0.8] 


0.142 


±0.033 


+0.057 
-0.051 



Table 3: Normalised differential cross section as a function of ^. 
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1 dE 

Mxd(i)* 




C>0.4 


i < 0.4, xi > 0.9 


C < 0.4, xi < 0.8 


i < 0.4, 0.8 <xi< 0.9 


(j)* (rad.) 


1 dE 
Mx d4>* 


(stat. ^ syst.) 


1 dE 

Mx d<p* 


(stat. syst.) 


1 dE 
Mxd<t>* 


(stat. syst.) 


1 dE 

Mx dtp* 


(stat. syst.) 


-3.016 


0.061 


+0.019 
-0.008 


0.122 


+0.014 
-0.017 


0.024 


+0.006 
-0.005 


0.072 


+0.008 

-0.012 


-2.765 


0.049 


+0.009 
-0.011 


0.288 


+0.026 
-0.024 


0.040 


+0.007 
-0.010 


0.114 


+0.015 
-0.011 


-2.513 


0.075 


+0.013 
-0.012 


0.335 


+0.029 
-0.043 


0.129 


+0.019 
-0.028 


0.239 


+0.023 
-0.018 


-2.262 


0.130 


+0.014 
-0.020 


0.139 


+0.020 
-0.014 


0.282 


+0.048 
-0.033 


0.289 


+0.020 
-0.022 


-2.011 


0.168 


+0.028 
-0.015 


0.050 


+0.010 
-0.007 


0.356 


+0.034 
-0.049 


0.151 


+0.012 
-0.023 


-1.759 


0.182 


+0.016 
-0.036 


0.031 


+0.005 
-0.006 


0.179 


+0.033 
-0.023 


0.077 


+0.010 
-0.008 


-1.508 


0.100 


+0.016 
-0.014 


0.027 


+0.007 
-0.006 


0.068 


+0.020 
-0.013 


0.054 


+0.008 
-0.009 


-1.257 


0.063 


+0.013 
-0.012 


0.021 


+0.005 
-0.007 


0.043 


+0.011 
-0.009 


0.040 


+0.008 
-0.007 


-1.005 


0.042 


+0.008 
-0.007 


0.026 


+0.006 
-0.006 


0.039 


+0.010 
-0.009 


0.029 


+0.011 
-0.004 


-0.754 


0.066 


+0.011 
-0.010 


0.043 


+0.007 
-0.007 


0.069 


+0.012 
-0.011 


0.072 


+0.009 
-0.011 


-0.503 


0.124 


+0.016 
-0.014 


0.089 


+0.019 
-0.012 


0.130 


+0.024 
-0.025 


0.118 


+0.013 
-0.018 


-0.251 


0.341 


+0.030 
-0.052 


0.325 


+0.030 
-0.040 


0.270 


+0.037 
-0.027 


0.276 


+0.034 
-0.024 


0.000 


0.661 


+0.060 
-0.050 


0.819 


+0.046 
-0.050 


0.512 


+0.046 
-0.046 


0.641 


+0.032 
-0.030 


0.251 


0.256 


+0.032 
-0.025 


0.303 


+0.048 
-0.027 


0.284 


+0.033 
-0.047 


0.333 


+0.024 
-0.034 


0.503 


0.106 


+0.019 
-0.014 


0.103 


+0.016 
-0.022 


0.124 


+0.023 
-0.024 


0.121 


+0.012 
-0.014 


0.754 


0.065 


+0.012 
-0.012 


0.041 


+0.008 
-0.008 


0.062 


+0.011 
-0.011 


0.045 


+0.009 
-0.006 


1.005 


0.035 


+0.006 
-0.008 


0.021 


+0.006 
-0.004 


0.035 


+0.010 
-0.006 


0.026 


+0.008 
-0.004 


1.257 


0.019 


+0.006 
-0.006 


0.025 


+0.006 
-0.007 


0.023 


+0.005 
-0.006 


0.036 


+0.005 
-0.007 


1.508 


0.022 


+0.005 
-0.007 


0.021 


+0.007 
-0.004 


0.029 


+0.007 
-0.006 


0.027 


+0.007 
-0.004 


1.759 


0.024 


+0.007 
-0.005 


0.037 


+0.006 
-0.007 


0.066 


+0.016 
-0.018 


0.042 


+0.008 
-0.006 


2.011 


0.039 


+0.011 
-0.007 


0.046 


+0.010 
-0.012 


0.233 


+0.034 
-0.055 


0.095 


+0.009 
-0.015 


2.262 


0.080 


+0.013 
-0.013 


0.106 


+0.012 
-0.026 


0.425 


+0.076 
-0.049 


0.195 


+0.016 
-0.026 


2.513 


0.352 


+0.031 
-0.048 


0.258 


+0.028 
-0.035 


0.363 


+0.064 
-0.039 


0.440 


+0.042 
-0.027 


2.765 


0.675 


+0.046 
-0.041 


0.439 


+0.058 
-0.026 


0.144 


+0.020 
-0.020 


0.335 


+0.021 
-0.033 


3.016 


0.249 


+0.027 
-0.027 


0.248 


+0.022 
-0.029 


0.051 


+0.010 
-0.012 


0.117 


+0.017 
-0.010 



Table 4: Energy-flow distribution in the three-jet plane, normalised to the total invariant mass, 
as a function of the azimuthal angle, 0*, in the different regions of the {C,,Xi) plane. 
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4 (pb) 


dr]*et 


^l 'p^' 


(stat.) (pb) 


(stat. syst.) (pb) 


[-3.6,-3.0] 


0.650 


±0.121 


+0.250 
-0.185 


[-3.0,-2.4] 


2.167 


±0.228 


+0.335 
-0.352 


[-2.4,-1.8] 


4.279 


±0.333 


+0.972 
-0.557 


[-1.8,-1.2] 


8.193 


±0.465 


+1.493 
-1.034 


[-1.2,-0.6] 


10.322 


±0.526 


+1.162 
-0.841 


[-0.6,0.0] 


9.384 


±0.506 


+0.940 
-0.958 


[ 0.0,0.6] 


7.132 


±0.429 


+0.704 
-0.856 


[ 0.6,1.2] 


8.359 


±0.445 


+1.295 
-0.837 


[ 1.2,1.8] 


9.417 


±0.471 


+1.171 
-0.857 


[ 1.8,2.4] 


9.543 


±0.605 


+1.241 
-1.338 


[ 2.4,3.0] 


1.039 


±0.145 


+0.381 
-0.203 



Table 5: Differential cross section for three-jet production as a function of the jet pseudorapidity, 
da/drj^^^, of the three jets. 



, -^^ (pb/GeV) 


dP|,j,t (GeV) 


4:. (pb/G^^) 


(stat.) (pb/GeV) 


(stat. syst.) (pb/GeV) 


[1.0,3.0] 


3.970 


±0.182 


+0.521 
-0.333 


[3.0,5.0] 


2.291 


±0.136 


+0.378 
-0.244 


[5.0,7.0] 


0.818 


±0.090 


+0.143 
-0.175 


[7.0,9.0] 


0.172 


±0.039 


+0.058 
-0.070 



Table 6: Differential cross section for three-jet production as a function of the jet transverse 



momentum, -Pxjet) of ^^e most-forward jet. 
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pM 




p{ip) (most-forward jet) 


p{f) (most-backward jet) 


dip (radians) 


Pi^) 


(stat.) 


(stat. syst.) 


Pif) 


(stat.) 


(stat. syst.) 


[0.0,0.2] 


1.524 


±0.057 


+0.086 
-0.083 


2.188 


±0.063 


+0.084 
-0.119 


[0.2,0.4] 


1.500 


±0.042 


+0.057 
-0.061 


1.273 


±0.038 


+0.127 
-0.054 


[0.4,0.6] 


0.877 


±0.029 


+0.043 
-0.053 


0.721 


±0.030 


+0.044 
-0.063 


[0.6,0.8] 


0.473 


±0.020 


+0.053 
-0.026 


0.350 


±0.017 


+0.024 
-0.034 


[0.8,1.0] 


0.300 


±0.015 


+0.022 
-0.019 


0.219 


±0.013 


+0.019 
-0.023 


[1.0,1.2] 


0.165 


±0.010 


+0.018 
-0.013 


0.125 


±0.009 


+0.015 
-0.011 


[1.2,1.4] 


0.103 


±0.008 


+0.010 
-0.015 


0.069 


±0.006 


+0.010 
-0.007 


[1.4,1.6] 


0.047 


±0.005 


+0.007 
-0.008 


0.025 


±0.003 


+0.008 
-0.005 



Table 7: The differential jet shape, p{ip), for the most-forward and most-backward jets with 
E*^^ > 9 GeV in three-jet events, where the Pomeron defines the forward direction. 
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